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Abstract —Social Spider Algorithm (SSA) is a recently pro¬ 
posed general-purpose real-parameter metaheuristic designed 
to solve global numerical optimization problems. This work 
systematically benchmarks SSA on a suite of 11 functions with 
different control parameters. We conduct parameter sensitivity 
analysis of SSA using advanced non-parametric statistical tests to 
generate statistically significant conclusion on the best performing 
parameter settings. The conclusion can be adopted in future work 
to reduce the effort in parameter tuning. In addition, we perform 
a success rate test to reveal the impact of the control parameters 
on the convergence speed of the algorithm. 

Index Terms —Social spider algorithm, global optimization, 
parameter sensitivity analysis, evolutionary computation, meta- 
heuristic. 

I. Introduction 

Social Spider Algorithm (SSA) is a recently proposed meta¬ 
heuristic designed to solve real-parameter black-box opti¬ 
mization problems m. SSA mimics the foraging behavior 
of the social spiders El for optimization. It has been ap¬ 
plied in solving both benchmark optimization problems m 
and real-world problems Em and demonstrated satisfactory 
performance. However, a parameter sensitivity analysis on the 
control parameters of SSA is still lacking. Although the op¬ 
timization performance of SSA in existing work outperforms 
the compared algorithms, it is significantly influenced by the 
four control parameters, i.e., population size, attenuation rate, 
change rate, and mask rate. Thus tuning these parameters 
properly is essential to demonstrate the efficacy of SSA as 
a problem solver. 

In general, the success of evolutionary algorithms (EAs) in 
solving optimization problems relies on the proper selection 
of control parameters 0- Control parameter tuning is the pro¬ 
cess where different combinations of parameters are assessed 
to achieve a best performance in terms of solution quality 
and/or computational complexity. Parameter sensitivity anal¬ 
ysis methods and applications have attracted much research 
attention in the past decades and many work has been done 
to fully exploit the searching power of EAs 0. An early 
example of such applications was conducted by Grefenstette 
Q, in which the control parameters of a genetic algorithm 
was tuned by another one so as to improve the performance 
on an image registration problem. Ever since then, parameter 
sensitivity tests have been performed on most of the existing 
EAs, including but not limited to PSO 0 and DE ||9l- 


In this work, we systematically study the impact of different 
parameter settings on a benchmark suite with a collection of 
optimization functions with different properties. We employ 
a set of 11 benchmark functions selected from a widely 
used benchmark function suite Qol and from the recent CEC 
single-objective continuous optimization testbed ifTTl . All test 
functions are the base functions from these two sources. We 
evaluate SSA using a wide range of parameter settings on the 
11 benchmark functions and present the best achieved results 
among all possible parameter combinations. In addition, we 
conduct a parameter sensitivity analysis with respect to all 
the combinations. The Eriedman test and the Hochberg post- 
hoc procedure inaiia are employed to generate statistically 
significant conclusion. The conclusion can be adopted in future 
related work to reduce the parameter tuning effort. 

The rest of this paper is organized as follows. Section II 
gives a brief review of SSA. Section III presents the experi¬ 
mental settings and simulation results. Section IV concludes 
this paper. 

II. Social Spider Algorithm 

A real-parameter black-box optimization problem aims to 
find the optimal decision variables, or the global best solution, 
to minimize/maximize one or multiple objective functions. In 
the optimization process, no knowledge of the function charac¬ 
teristics is known to the problem solver. SSA is a population- 
based general-purpose stochastic optimization problem solver 
and has been shown to be effective and efficient in solving 
real-parameter black-box optimization problems m . It has also 
been applied in scientific and engineering research 00. 

In SSA, the search space of an optimization problem is 
considered as a hyper-dimensional spider web, on which a 
number of artificial spiders can move freely. Each position on 
the web corresponds to a feasible solution to the optimization 
problem. In addition, the web also serves as a transmission 
media for the vibrations generated by the movement of the 
spiders. Details of SSA will be introduced in the following 
subsections. 

A. Spider 

The artificial spiders are the basic operating agents of 
SSA. Each spider possesses a position on the hyper-dimension 
spider web, and the fitness value of this position is assigned 
to the spider. Each spider holds a memory storing its status as 
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well as optimization parameters, namely, its current position, 
current fitness value, following vibration at previous iteration, 
inactive degree, previous movements, and dimension mask. All 
these pieces of information are utilized to guide the spider to 
search for the global optimum, and will be introduced later. 

B. Vibration 

In SSA, a spider will generate a vibration upon moving to 
a new position different from the one in the previous iteration. 
Vibration is a unique feature of SSA, which distinguishes SSA 
from other meta-heuristics. The vibrations serve as medium 
for the spiders to perform lossy communications with others 
to form a collective social knowledge. 

A vibration in SSA is characterized by two properties, 
namely, its source position and source intensity. When a spider 
moves to a new position P, it generates a vibration whose 
source position is P. The source intensity of this vibration is 
calculated using the following equation: 

where I is the source intensity, f{P) is the fitness value of P, 
and C is a constant smaller than the infimum of the objective 
functiorQ 

After a vibration is generated, it will propagate over the 
spider web. As a form of energy, vibration attenuates over 
distance. This physical phenomenon is considered in SSA. 
The vibration attenuation process is defined with the following 
equation: 

I‘^ = Ixexp {-^—), ( 2 ) 

a XTa 

where is the attenuated intensity after propagating for 
distance d, a is the average of the standard deviation of 
all spider positions along each dimension. In this equation, 
Ta S (0, +oo) is introduced to control the attenuation rate of 
the vibration intensity, and will be tested as a user-controlled 
parameter in this paper. 

C. SSA Search Pattern 

SSA manipulates a population of artificial spiders via a 
series of optimization step. Specifically, each iteration of SSA 
can be divided into the following steps. 

1) Fitness Evaluation: At the beginning of each iteration, 
the fitness values of the positions possessed of all spiders in 
the population is re-evaluated. These fitness values are later 
utilized in the vibration generation and propagation process. 

2) Vibration Generation: After all fitness values are calcu¬ 
lated, each spider in the population will generate a new vibra¬ 
tion at its current position. Its source intensity is calculated 
using Q. This vibration is immediately propagated to all the 
spiders in the population via the spider web with attenuation. 
Thus, each spider will receive \pop\ different vibrations from 
the web when all spiders have generated new vibrations in the 

*We consider minimization problem in this paper. For maximization prob¬ 
lems, C is a constant larger than the supremum of the objective function, and 
the denominator of the equation is C — f{P)- 


current iteration, where \pop\ is the population size. Upon the 
receipt of these \pop\ vibrations, the spider will select the one 
with the largest attenuated vibration intensity, and compare 
it with the one it stored, called the “following vibration” in 
the previous iteration. The vibration with a larger intensity 
is stored as the updated following vibration in the current 
iteration. If the spider does not change its following vibration, 
its inactive degree is incremented by one. Otherwise this 
degree is reset to zero. This ends the vibration generation step. 

3) Mask Changing: When all spiders have determined their 
following vibration, the algorithm proceeds to manipulate their 
positions on the spider web. Here SSA utilizes a dimension 
mask to introduce randomness in the following movement 
process. Each spider holds a dimension mask which is a 
0-1 binary vector of length D, where D is the number of 
dimensions of the optimization problem. In each iteration, the 
spider has a probability of 1 — pf" to change its mask where 
Pc € (0,1) is a user-controlled parameter this paper aims to 
study. Here we use dm to denote the inactive degree of the 
spider. If the mask is decided to be changed, each bit of the 
mask has a probability of Pm to be assigned with a one, and 
1 — Pm to be a zero, pm is another user-controlled parameter 
in the parameter sensitivity analysis. In this process, each bit 
is changed independently, except for the case where all bits 
are set to zero. In such a case, one random bit of the mask is 
set to one in order to avoid premature convergence. 

4) Random Walk: With the following vibration and dimen¬ 
sion mask determined, each spider then generates a following 
position based on these two pieces of information and all re¬ 
ceived vibrations. This following position is employed together 
with the spider’s current position to generate a new position on 
the spider web. This new position is assigned to the spider and 
this ends the random walk step. A more detailed formulation 
can be found in m. 

D. Control Parameters 

In SSA, four user-controlled parameters shall be fine-tuned 
in order to achieve good performance: 

1) Population Size \pop\: The size of the spider population 
determines the individual diversity and influences the conver¬ 
gence speed. The selection of this parameter greatly influences 
the computational complexity as the standard deviation of 
all solutions shall be calculated when generating attenuated 
vibration intensities using In addition, while a larger 
population may increase the diversity of the population which 
in turn increases the exploration ability, a smaller population 
may result in better convergence speed. 

2) Attenuation Rate ra-' This parameter controls the attenu¬ 
ation strength imposed on the vibration intensities. A larger 
will make the exponential part of Q approach one, resulting 
in a weaker attenuation on the vibration intensity. This means 
that the communication is less lossy and best positions in the 
population are more likely to be broadcast to all spiders. On 
the contrary, a smaller Va results in a stronger attenuation, 
rendering the information propagation more lossy. This may 




TABLE I 

Benchmark Functions 


Function 

Search Range 

Dimension (D) 

Name 

/iW = 

i = l 

[-100,1001^ - 0 

30 

Generalized Sphere Function 

/ 2 (x) = 

i=2 

[-100,1001^ - 0 

30 

Generalized Cigar Function 

D i 

/3(x) = 

[-100,1001^ - 0 

30 

Schwefel’s Function 1.2 

i = l j = l 




fi{x) = maxjxil 
i 

[-100,1001^ - 0 

30 

SchwefeTs Function 2.21 


D-l 


/sW = - 1)^) 

i = l 

D 

/g(x) =418.9829D -'^g(xi + 420.9687) 

i=l 

{ Zisin{\zi\^/‘^) if|2i|<500 

Vi sin(\/Ii/d) - 

Vi sin(x/Ij/^) - if < -500 


[—100,100]^ — o 30 Generalized Rosenbrock’s Function 


[—1000,1000)^ — lOo 30 Modified SchwefeTs Function 


Vi = 

D 


500 - mod(zi, 500) if Zi > 500 

mod(- 2 i, 500) - 500 if Zi < -500 


frix) = — 10cos(27r3ij) + 10) 

i = l 


[-5.12,5.12)° - 0.0512O 

30 

Generalized Rastrigin’s Function 

/8(x) = -20exp(-0.2^ 

^ n 1 ^ 

— xf) — exp[— cos(27rx 

n •“ n •“ 

1=1 1=1 

)] 

[-32,32)° - 0.32O 

30 

Ackley’s Function 

-l-20-l-e 






/g(x) = - xf — TT cos(^) -F 1 

^ ^ 4000 ^ * y wi 


[-600,600)° - 6o 

30 

Generalized Griewank’s Function 

/io(x) =g{xi,X2) +g{x2,X3) H- \- g(xn--i,Xn) + g{x 

n,Xl] 




g{x,y) = 0.5 + 

10 ^ 

i=l j 

(sin^(x/a;3 + y'^) — 0.5) 


[-100,100)° - 0 

30 

Scaffer’s Function F6 

(l-|-0.001(x2 +i/2))2 
^j2^Xi -round{2Jxi)| UTTt 

10 

[-5.12,5.12)° - 0.0512O 

30 

Katsuura’s Function 

-J “ 2-? 

£>2 


potentially obstruct the algorithm from fast convergence but 
prevent it from premature convergence. 

3) Mask Changing Probabilities pc and pm-' These two 
parameters cooperate to determine the dimension mask of each 
spider in the population. A large pc will result in frequent 
mask changing, which can be beneficial to solve optimization 
problems with a great number of local optima. A large pm 
will result in more I’s in the dimension mask, which means 
that more noise is introduced in the following position when 
performing random walk ||T|. 

III. Experimental Settings and Simulation Results 

We evaluate the optimization performance of SSA with 
all possible 900 combinations of six \pop\ values, six 
values, hve pc values, and five pm values on 11 numerical test 
functions selected from the base functions of a widely-adopted 
benchmark suite Go) and a recent CEC single-objective 
continuous optimization testbed GD- These functions are 
employed to construct the competition functions in recent 


CEC single-objective optimization competitions GUGIGS, 
and are considered representative for covering most possible 
features of real-world application problems. In addition, we re¬ 
port the simulation results of the best performing combination 
of control parameters with additional performance analysis. 
To provide a statistically signihcant and reliable guideline 
for parameter selection, we perform the parameter sensitivity 
analysis using a rigid non-parametric statistical test framework 
to compare all 900 test cases. 

A. Benchmark Functions 

In this work, we adopt 11 base functions from Ga and 
from the CEC 2013 single-objective continuous optimization 
testbeds. These base functions are grouped into two classes: 
uni-modal functions / 1-/4 and multi-modal functions / 5 -/ 11 . 
The benchmark functions are listed in Table U 

As one of the objective is to test the convergence ability of 
the algorithm, we use a novel method to shift the position of 
the global optima yet maintain the accuracy of computation to 













smaller than 10 m. Previously, most benchmark function 
suites move the global optima from 0 to a random position 
in order to avoid unfair competition by algorithms that only 
exploits the center of the search space. However, due to the 
storage method in modem computer, the accuracy of a floating 
point number is limited to its size, which makes it difficult to 
store large numbers with small tails, e.g. a = 50 + 10“®°. In 
order to preserve a high accuracy and avoid global optima at 
0 , the benchmark functions in this paper adopt a search space 
shift vector o on the search spaces. For example, /i originally 
is designed to search for global optimum in [—100,100]^. 
In our benchmark function suite, the search space for /i is 
modified to [—100,100]^ — o = [—100 — oi, 100 — oi] x 
[—100 — 02 ,100 — 02 ] X • • • X [—100 — 0£), 100 — od]- By this 
means, the global optimum, although still at 0 , is not the center 
of the search space anymore. The value of o is set according 
to Cl. 

B. Experimental Settings 

It is a very time consuming task to choose a proper parame¬ 
ter configuration for SSA to generate outstanding performance. 
As many real-world application problems are computationally 
expensive, it may take much time to finish one run of the 
algorithm, which makes the trial-and-eiTor method for param¬ 
eter tuning impractical. One solution to this problem is to 
employ the “fixed scheme” for parameter configuration, which 
uses a fixed parameter setting throughout the search based on 
extensive empirical studies on a wide range of test functions 
with diverse properties. Such empirical studies can deduce 
guidelines for choosing parameter configurations leading to 
consistently satisfactory performance over most optimization 
problems. 

In this paper, we employ advanced non-parametric statistical 
tests to conduct parameter sensitivity analysis and to determine 
some parameter combinations that can lead to statistically sig¬ 
nificant better performance than others. Our statistical analysis 
is a two-step framework. In the first step, the Friedman test 
MM is conducted to compare 900 parameter settings over 
the 11 test functions to test the null hypothesis that no test 
case performs significantly different from others. If this null 
hypothesis is rejected, a Hochberg post-hoc procedure ifTBlUfTTll 
is adopted to determine which settings perform significantly 
better than others. 

We test SSA with 900 parameter combination settings, 
namely [\pop\,ra,Pc,PTn] G [10,20,30,40,50,70] x [0.2, 
0.5,1.0,2.0,4.0,8.0] X [0.1,0.3,0.5,0.7,0.9] x [0.1,0.3,0.5, 
0.7,0.9] on all 11 benchmark functions on 30 dimensions, 
respectively. For each benchmark function, SSA with each of 
900 parameter settings are executed 51 times. The stopping 
criteria for all runs are set to 300 000 function evaluations, 
according to the suggestion by M- 

C. Parameter Sensitivity Analysis 

Table [II] reports the results of the parameter sensitivity anal¬ 
ysis. In the table, different values are presented horizontally 
while Pc X Pm combinations are presented vertically. If a 


TABLE II 

Parameter Sensitivity Analysis Result 




Va 

Pm 

Pc 

0.2 

0.5 

1.0 

2.0 

4.0 

8.0 


0.1 

20/30 

20/30 

20/40 

20/30 


- 


0.3 

20/30 

20/30 

20/30/40 

20/30 


- 

0.1 

0.5 

10/20/30 

20/30 

20/30/40 

20/30/40 


- 


0.7 

10/20/30 

10/20/30 

20/30/40 

20/30/40 


- 


0.9 

10/20/30 

10/20/30 

10/20/30 

20/30 

- 

- 


0.1 






- 


0.3 






- 

0.3 

0.5 

20/30 

20/30 

20 

20/30/40/50 


- 


0.7 




40 


- 


0.9 

- 

- 

- 

- 

- 

- 


0.1 






- 


0.3 






- 

0.5 

0.5 





30 

- 


0.7 






- 


0.9 

- 

- 

- 

- 

- 

- 


0.1 






- 


0.3 






- 

0.7 

0.5 






- 


0.7 






- 


0.9 

- 

- 

- 

- 

- 

- 


0.1 



70 



- 


0.3 






- 

0.9 

0.5 






- 


0.7 


40/50/70 




- 


0.9 

- 

50/70 

70 

- 

- 

- 


specific \pop\ x Va x pc x pm combination performs signifi¬ 
cantly better than others, \pop\ is labeled in the corresponding 
'fa>^Pc>^ Pm cell, i.e., all parameter combination labeled in 
the table performs significantly better than those not labeled. 

From this table, we have the following observations: 

• Medium \pop\ values (e.g., 20, 30, and 40), small to 
medium ra values (e.g., 0.5, 1, and 2), medium pc values 
(e.g., 0.5 and 0.7), and small pm values (e.g., 0.1 and 
0.3) are preferred by the statistical analysis. 

• A roughly most preferred parameter combination for 
solving the benchmark functions in this paper is 
[\pop\,ra,Pc,Pm] = [30/1.0/0.7/0.1]. This result ac¬ 
cords with our preliminary observation in IT]. 

• There are some noisy data points at pm = 0.9. After 
a further investigation of the raw simulation data, we 
found that the parameter combinations with pm = 0.9 
sacrifices the performance on uni-modal functions for 
some tiny improvements in /n, which is multi-modal, 
non-separable, and the number of local optima is huge. 
As the improvement is barely enough to outperform 
others, it possesses advantage in the statistical analysis. 
However, they are not preferred when considering the 
overall performance on all benchmark functions. 


The simulation results for the recommended parameter 


combination by the statistical test are listed in Table III 


Besides the numerical analysis on the simulation results, 
convergence speed is also a critical performance metric for 
meta-heuristics. In order to investigate the impact of parameter 
settings on the convergence speed, we adopt the success rate 
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Fig. 1. Success Rate Plots of Different Parameters. 


analysis ||9l to see the influence of the four parameters. In this 
success rate analysis, a simulation run is considered successful 
when the best fitness value is smaller than 10“®. For each run, 
the best fitness value is recorded every 3000 FEs and there are 
in total 100 “success bits” per run, recorded whether the run 
is successful up to the specific FEs. Thus 576 (11 x 51) runs 
of success bits for each parameter setting can be accumulated 
and presented as a plot of success rate over all test functions. 

In this test we employ the previously recommended param¬ 
eter setting [30/1.0/0.7/0.1] as the standard setting, denoted 
as “Standard”. Different values of the four parameters are 
employed to compare with the standard setting, and only one 
parameter is consider in each comparison. The success rate 
plots are shown in Eigure In this figure, plot (a) studies 
the impact of \pop\ values on the convergence speed, plot (b) 
studies the impact of r^, plot (c) studies the impact of pc, and 
plot (d) studies the impact of pm- Erom this figure, we have 
the following observations; 

• The impact of \pop\ on the convergence speed is roughly 
the most significant one. This observations accords with 
our previous analysis in Section II-D. 

• Reasonable r^, Pm, and pc values yield similar conver¬ 


gence speed. Eor instance, the convergence plots for ra € 
[1.0, 2.0] are almost identical for the whole process, and 
the plot for = 0.2 joins them at approximately 35% 
of the whole optimization process. However, = 8.0, 
which is a very rare setting for this parameter, cannot 
reach the success rates of others for the whole process. 

IV. Conclusion 

In this work we systematically evaluate the performance 
of SSA on a benchmark suite of 11 different func¬ 
tions. We carry out a parameter sensitivity analysis on 
900 possible parameter settings on the four control pa¬ 
rameters ([\pop\,ra,Pc,PTn] G [10,20,30,40,50,70] x [0.2, 
0.5,1.0,2.0,4.0,8.0] X [0.1,0.3,0.5,0.7,0.9] x [0.1,0.3,0.5, 
0.7, 0.9]) using a rigid statistical test framework. The results 
indicate that medium \pop\ values (e.g., 20, 30, and 40), small 
to medium Va values (e.g., 0.5, 1, and 2), medium pc values 
(e.g., 0.5 and 0.7), and small pm values (e.g., 0.1 and 0.3) can 
lead to statistically significantly better performance than other 
settings. This is a reliable rule of thumb for parameter tuning 
of SSA in the future work. In addition, we using the success 
rate test to analyze the impact of these four parameters on the 











































TABLE III 

Parameter Sensitivity Analysis Result 


Function 

Mean 

Std. Div. 

Best 

Worst 

Median 

h 

1.1321E-73 

1.3464E-73 

2.4108E-75 

7.8390E-73 

7.7373E-74 

h 

9.2887E-70 

5.2257E-69 

3.3337E-78 

3.7499E-68 

4.3240E-74 

h 

1.3651E+00 

8.9476E-01 

1.5132E-01 

3.8509E+00 

1.1759E+00 

U 

3.9103E-03 

4.4044E-03 

2.7803E-04 

2.523 lE-02 

2.1600E-03 

/s 

6.2745E+01 

1.0091E+02 

4.0258E-04 

5.1402E+02 

1.8734E+01 

/6 

1.9676E+01 

1.1145E+02 

1.6694E-01 

8.0221E+02 

1.4719E+00 

/t 

6.1569E-07 

4.1956E-06 

O.OOOOE+00 

3.0274E-05 

O.OOOOE+00 

fs 

7.4798E-15 

4.9258E-16 

3.9968E-15 

7.5495E-15 

7.5495E-15 

h 

6.2035E-13 

1.8712E-12 

O.OOOOE+00 

1.2479E-11 

3.9191E-14 

fio 

6.4117E+00 

3.9909E-01 

5.2466E+00 

7.1609E+00 

6.4479E+00 

hi 

2.1116E-01 

1.361 lE-01 

2.5308E-03 

4.5302E-01 

2.2008E-01 


convergence speed of SSA. The simulation results indicate that 
while \pop\ value has a significant influence, rational pm, 
and Pc values generally lead to similar convergence speed. 
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